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摘要 :目的 研发 一 种 基于 高 维 全 变 分 (high dimension total variation, HDTV) 先 验 的 CT 心肌 灌注 (computed tomography 


myocardial perfusion, CT-MP) 去 卷 积 方法 。 方 法 


先 对 低 剂量 CT-MP 数 据 进行 灌注 去 卷 积 建 模 ,其 后 通过 引入 HDTV 正则 化 
修正 模型 最 小 化 解 的 一 致 性 。 其 中 ,HDTV 利用 心肌 血管 空间 结构 相似 性 和 心肌 血 流 信号 时 间 连 续 性 。 结 果 XCAT 心 肌 灌注 


仿真 数据 和 猪 心 肌 灌注 数据 实验 表明 , 同 已 有 的 方法 相 比 ,本 方法 能 更 有 效 地 抑制 低 剂量 血 流动 力学 参数 图 中 的 噪声 和 伪 影 ， 
同时 较 好 地 保持 具有 诊断 意义 的 结构 信息 。 结 论 本 文 方法 可 实现 低 剂量 CT-MPI 血 流动 力学 参数 图 的 优质 成 像 。 
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Robust low-dose CT myocardial perfusion deconvolution via high-dimension total 


variation regularization 
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Abstract: Objective To develop a computed tomography myocardial perfusion (CT-MP) deconvolution algorithm by 


incorporating high-dimension total variation (HDTV) regularization. Methods A perfusion deconvolution model was 
formulated for the low-dose CT-MPI data, followed by HDTV regularization to regularize the consistency of the solution by 


fusing the spatial correlation of the vascular structure and the temporal continuation of the blood flow signal. Results Both 
qualitative and quantitative studies were conducted using XCAT and pig myocardial perfusion data to evaluate the present 
algorithm. The experimental results showed that this algorithm achieved hemodynamic parameter maps with better 


performances than the existing methods in terms of streak-artifacts suppression, noise-resolution tradeoff, and diagnosis 
structure preservation. Conclusion The proposed algorithm can achieve high-quality hemodynamic parameter maps in 


low-dose CT-MPI. 
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CT 心肌 灌注 成 像 (CEMPD 作 为 一 种 有 效 的 非 侵 
人 功能 成 像 方法 ,其 在 冠状 动脉 疾病 (CAD ) 的 诊断 和 危 
伶 分 级 方面 有 显著 的 优势 "。 在 CTMPI 的 研究 中 ,可 
以 通过 心肌 灌注 数据 估计 血 流动 力学 参数 图 ,比如 心肌 
Hii (MBF) ,心肌 血 容量 (MBV) 和 平均 通过 时 间 
(MIT) ,这 些 血 流 动力 学 参数 图 可 以 为 临床 医生 诊断 
提供 重要 辅助 信息 。CT 血管 造影 技术 和 CTMPI 技 术 
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联合 使 用 ,不 仅 能 提供 冠状 动脉 狭 罕 的 功能 信息 而 且 还 
能 获得 其 狭窄 部 位 的 解剖 结构 信息 ”。 根 据 扫描 协议 ， 
CT-MPI 对 感 兴趣 层面 连续 多 次 增强 扫描 ,和 传统 的 CT 
扫描 相 比 ,这 大 大 增加 了 病人 所 受 的 罩 射 剂量 。 大 量 研 
究 表明 ,过 高 的 X 射 线 照 射 会 诱发 癌症 ,白血病 或 其 他 
遗传 性 疾病 ”。 当 前 ,临床 上 采用 的 低 剂 量 扫描 方案 虽 
然 是 一 种 降低 辐射 剂量 简单 易 行 的 手段 ,但 常常 会 导致 
所 采集 的 投影 数据 受到 量子 噪声 的 影响 ,使 得 重建 CT 
心肌 灌注 图 像 质量 出 现 严重 退化 ,从 而 影响 诊断 的 精确 
度 。 因 此 ,在 保证 一 定 图 像 质量 的 条 件 下 ,如 何 降低 
CTMPI 中 X 射 线 辐射 剂量 已 成 为 当前 CT 成 像 尤 其 是 
CT-MPI 领 域 最 为 重要 的 课题 之 一 。 

当前 ,针对 CTMPI 剂 量 降低 问题 ,除了 采用 心 电 
门 控 技术 和 优化 扫描 模式 ,许多 抑制 图 像 噪声 的 方法 相 
继 被 提出 *。 如 Speidel 等 在 不 考虑 运动 伪 影 的 条 件 
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下 采用 HYPR 技 术 恢 复 低 剂 量 心肌 灌注 图 像 ,该 技术 能 
在 一 定 程度 上 抑制 CT 心肌 灌注 图 像 噪声 。Ritman 等 
考虑 了 由 心脏 运动 带 来 的 运动 伪 影 影响 ,利用 多 排 探 
测 右 CT 进行 部 分 扫描 重建 的 MPI 人 研究 ,提出 了 一 种 基 
于 目标 的 空间 频率 滤波 策略 ,该 策略 能 够 减少 因 部 分 
扫描 引入 的 重建 伪 影 。Tao 等 "将 统计 迭代 (statistical 
iterative reconstruction, SIR) 方 法 应 用 于 CT-MPI, 这 种 
直接 作用 于 投影 数据 的 CT 迭代 重建 方法 能 有 效 的 抑制 
噪声 和 条 行 伪 影 ;Bian 等 "提出 了 基于 SIR 的 运动 自 适 
应 稀 玻 先 验算 法 应 用 于 低 剂量 CT 心肌 灌注 图 像 重建 ， 
该 方法 利用 了 CT 心肌 灌注 图 像 的 空间 及 时 间 结 构 的 稀 
玻 特 性 并 将 运动 校正 过 程 引 入 到 迭代 重建 过 程 中 ,该 方 
法 在 噪声 消除 , 伪 影 抑制 和 细节 保持 方面 均 有 恨 好 效 
果 。 昌 然 上 述 方法 能 够 获得 优质 的 CT 心肌 灌注 图 像 ， 
但 是 这 些 方法 在 利用 得 到 的 CT 心肌 灌注 网 像 估计 血 流 
动力 学 参数 图 时 采用 不 稳定 的 截断 奇异 值 分 解 去 卷 积 
方法 。 根 据 文献 [9, 10] , 若 要 获得 具有 诊断 意义 的 血 流 
动力 学 参数 图 ,必须 采用 和 鲁 棒 的 灌注 去 卷 积 方法 。 
此 ,为 了 提高 血 流动 力学 参数 图 的 估计 精度 ,近年 来 许 
多 灌注 去 卷 积 方法 相继 被 提出 ,Frindel 等 "构建 了 
基于 时 空 正 则 化 去 卷 积 模型 应 用 于 脑 灌注 加 权 MRI 
(perfusion-weighted magnetic resonance imaging, 
PW-MRI) ,并 提出 了 全 局 收敛 算法 求解 相应 的 目标 函 
数 。Fang 等 "提出 基于 全 变 分 (TV) 正 则 化 去 卷 积 方法 
应 用 于 低 剂量 CT 脑 灌注 成 像 ,该 方法 能 够 有 效 改 善 脑 
血 流 动力 学 参数 图 。 

受 基 于 TV 正则 化 去 卷 积 模型 启发 ,本 文 提出 一 种 
基于 高 维 全 变 分 (High dimension total variation, 
HDTV) 正 则 化 的 CT 心肌 灌注 去 卷 积 方法 。 其 中 ， 
HDTV 正则 化 利用 了 心肌 血管 空间 结构 相似 性 和 心肌 
血 流 信号 时 间 连 续 性 。 同 时 ,针对 新 构建 的 CT-MPI 去 
卷 积 目标 函数 的 求解 ,本 文 提出 一 种 改进 的 交 蔡 方向 法 
优化 策略 。 心 肌 灌 注 仿真 数据 和 猪 心 肌 灌 注 数 据 实验 
表明 ,与 传统 的 去 卷 积 方法 "“ 相 比 ,本 文 方法 能 更 有 效 
地 抑制 低 剂 量 血 流动 力学 参数 图 中 的 噪声 和 伪 影 ,同时 
较 好 地 保持 具有 诊断 意义 的 结构 信息 。 


1 方法 
1.1 CT 灌注 参数 估计 的 去 卷 积 法 

基于 示 踪 剂 稀释 理论 ”和 CT 灌注 参数 的 计算 
模型 ,组 织 剩余 函数 R(1) 计 算 方式 如 下 : 

C(t) =MBF: C,(t)*R(t) =C,(t)*k(t) (1) 
其 中 ,x 为 卷 积 算 子 符号 ,Ci(1) 表 示 动 脉 对 比 剂 浓 度 ,也 
称 之 为 动脉 输入 函数 (arterial input function, AIF),C(7) 
表示 组 织 内 的 对 比 剂 浓度 ,K(D = MBF - R(t); (1) 式 中 
C(7) 和 Cl(7?) 可 从 灌注 数据 直接 测量 获得 ,k(1) 是 未 知 的 
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HEZE, MBF, MBV 和 MTT 等 灌注 参数 可 以 通过 K(1) 
求 得 。 根 据 文献 L10, 14] , 卷 积 方程 式 (1) 可 被 离散 化 为 
和 矩阵 形式 ,对 于 感 兴趣 层面 单个 像素 的 卷 积 离散 化 矩阵 
表达 式 为 : 

Ak=c (2) 

此 时 A eR”*”,N 是 采样 次 数 ,ceR" 和 keR"。 相 
对 于 (2) 式 ,整个 感 兴趣 层面 的 卷 积 离散 化 矩阵 表达 式 为 : 

AK=C (3) 
其 中 C= [crc] ER K= [krki] ERY”, LEIN 
兴趣 层面 的 像素 个 数 ,(1) 式 是 基于 这 样 的 假设 :对 比 剂 
FIA ATF 的 时 间 不 小 于 到 达 感 兴趣 层面 组 织 的 时 间 ; 实 
际 中 ,我 们 测 得 的 A 正 并 不 一 定 是 该 组 织 的 AIF ,一般 
情况 下 AIF 都 会 延迟 于 C(D),AIE 和 C() 之 间 的 时 间 漂 移 
会 导致 计算 的 灌注 参数 不 精确 ;为 了 避免 对 比 剂 延 
述 问 题 ,许多 研究 提出 了 不 同 的 解决 方法 ,其 中 
Wu 等 “采用 块 循环 矩阵 (4i,) 去 卷 积 法 效果 较为 显著 ， 
A 中 PTCA circ 表示 如 下 : 


(deirc) = | Atcan(tijs 小 forj<i a 


= 


4 


a Atan tusia) forj>i 
这 里 ,At 是 数据 采样 时 间 间 隔 ,i,j= 1,…N,M=2N， 
和 矩阵 E 和 C 也 要 从 N 维 扩展 到 M 维 ,将 (3) 式 替换 为 : 
HX=Y (5) 
ERTH =A oiro X E RFI © RR 分 别 是 对 K 和 C 补 
零 扩 展 后 的 矩阵 "9。 
1.2 基于 HDTV 正则 化 的 去 卷 积 模型 
关于 (5) 式 的 最 小 二 乘 解 等 价 于 极 小 化 该 式 的 二 次 
欧 几 里 得 残 差 范 数 ,可 计算 如 下 : 
Xs = arg min (|| HX- na) (6) 
KeRY*r 


其 中 , 互 是 一 个 病态 的 托 普 利 菊 矩阵 ,因此 和 并 不 是 一 
个 正定 的 最 小 二 乘 解 ,只 要 7 中 微小 的 变化 (投影 噪声 ) 
就 会 导致 如 很 大 变动 ;为 了 得 到 优质 的 X ,考虑 到 心肌 
血管 空间 结构 相似 性 和 心肌 血 流 信号 时 间 连 续 性 ,本 文 
提出 了 基于 HDTV 正则 化 的 去 卷 积 模型 ,可 表示 为 : 


ee arg min{ AHX- valk +iy(X)) (7) 
eR” 
FEP, EMNET (x yz LY 
2 
v(xX)=| Sor xet (8) 
i=0 


t 
HP D =D,,V%o=1,0,=D,,¥,=10,=1¥,=D,;D,, 
D, 和 D 是 分 别 沿 着 rz,y 和 (三维 的 有 限 差分 算 子 ; 为 正则 
化 参数 。 
1.3 优化 策略 

为 了 极 小 化 (7) 式 ,通过 引入 辅助 变量 $ 把 (7) 式 转 
变 成 一 个 等 价 的 约束 问题 ,如 下 : 
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seed 2 _ 
min || HX-Y|[, +4 i SS 


(9) 


t 


2 2 
>is 
i=0 

©; X¥;; i=0,1,2 


对 于 (9) 式 求解 等 价 于 极 小 化 下 面 的 增 广 拉 格 明 日 
+ 


me” 
2 
之 | sf 
i=0 (10) 


2 olor xe, + 1(Z,07X¥,-5) 
KE ZONA AEF BARETI BN, PE A 
由 定义 为 (4,B)=trace( A7B)。 我 们 使 用 交替 极 小 化 方 
法 极 小 化 (10) 式 ,具体 如 下 : 

Xr1= arg min|| HX- 7|| + PSyorxy | (11) 


1 2 
U(X,S)) =5||HX- Yl, +4 


forse, + 4).s|| 


= le 
Seo1= ang min > > 
(12) 
; 1=0,1,2 


i=0 
>I si 
i= 0 


A {Of IE-S i=0,1,2 (13) 


AS SCR FASEB ES EAE C11) SK, C12) SUR AE 
使 用 多 维 收缩 算法 


; A, 
图 1 XCAT 心 肌 灌注 图 像 和 低 剂 量 仿真 心肌 灌注 CT 图 像 
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P; 2 
Sini= 3 (Se é (14) 
i=0 


Sof ’ 


i= 


ae 


此 时 Pp, = OT XW, + A (13) 式 是 对 拉 格 朗 日 乘 子 2 
的 更 新 。 


2 结果 与 讨论 

为 了 验证 新 方法 在 低 剂量 CT 心肌 灌注 去 卷 积 
的 有 效 性 ,我 们 通过 XCAT 体 模仿 真 和 动物 数据 实验 分 
别 进行 了 定性 分 析 和 定量 分 析 。 
2.1 XCAT 心 肌 灌注 体 模 实验 

本 实验 中 ,选择 修订 的 XCAT™ 体 模仿 真 生 成 CT 
心肌 灌注 数据 。CT 仿 真 扫描 的 成 像 参 数 如 下 : XCAT 
体 模 图 像 大 小 为 512x512, 层 厚 为 1 mm; 射 源 到 旋转 中 
心 和 探测 器 的 距离 分 别 为 146 mm 和 1361 mm; 管 电压 
为 120 kVp, 旋 转角 在 [0,2x] 间 采样 值 为 1160, 探 测 器 个 
数 为 672。 低 剂量 投影 数据 生成 方法 参考 文献 [22-23 ] ， 
其 中 光子 总 数 设 为 2.5x105, 噪 声 方差 c2 = 10, 最 后 对 噪 
声 测量 数 据 进 行 EBP 重 建 获得 低 剂量 CTEMPI 图 像 ;为 
了 体现 本 文 方法 对 缺 血 心肌 显示 的 优越 性 ,我 们 通过 临 
床 医生 指导 在 XCAT 体 模 上 定义 了 一 块 马蹄 形 缺 血 心 
肌 区 域 ;图 1 给 出 了 真 值 和 低 剂 量 仿真 XCAT 心肌 灌注 
图 像 ,其 中 黑色 箭头 所 示 为 缺 血 心肌 。 


Fig.1 XCAT myocardial perfusion image and low-dose image reconstructed from the 


simulated projection data. A: XCAT myocardial perfusion phantom; B: low-dose 
myocardial perfusion image reconstructed from the simulated projection data. 


为 了 定量 分 析 不 同方 法 的 血 流动 力学 参数 图 像 质 
量 , 本 文 计算 了 标准 参数 图 像 与 四 种 去 卷 积 方法 得 到 的 
参数 图 像 之 间 的 峰值 信 噪 比 (peak signal to noise ratio, 
PSNR) 和 相对 均 方 根 误差 (relative root mean square 
error, IRMSE) ,其 计算 公式 如 下 : 


MAX? (Hxwue) 


PSNR=10 logio (15) 


Q 


Dhala) -p(n) (Q- 1) 


n=1 


Suln) - Hyteue(”)) 


rRMSE = |2=" (16) 


Q 2 
> (wln)) 


Ausaw 表 示 标 准 参数 图 像 ,表示 低 剂量 去 卷 积 得 到 
的 参数 图 像 ,0 是 感 兴趣 层面 内 的 像素 个 数 。 

图 2 所 示 为 不 同 去 卷 积 方法 得 到 XCAT 体 模 心肌 
灌注 MBF,MBV 和 MTT 血 流动 力学 参数 图 像 , 其 中 这 
三 种 去 卷 积 方法 :标准 截断 奇异 值 分 解 (standard 
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truncated singular value decomposition, sSVD) JNA"! 
快 循环 矩阵 截断 奇异 值 分 解 (block-circulant truncated 
singular value decomposition, bSVD) 方 法 2 和 Tikhonov 
正则 化 方法 “目前 广泛 应 用 于 商业 灌注 CT 设备 中 。 从 
图 2 可 以 看 出 ,由 sSSVD 方 法 得 到 的 三 种 参数 图 均 含 有 
严重 的 噪声 和 条 形 伪 影 ,正常 心肌 和 缺 血 心肌 区 域 几 乎 
被 噪声 淹没 ;bSVD Ail Tikhonov 的 结果 稍 有 改善 ,但 是 
与 Phantom 参数 图 相 比 存在 较 大 差距 , ACHE MCE 
MBF, MBV 和 MTT 参 数 图 像 中 分 辩 缺 血 心 肌 轮 廓 ;本 
文 方法 结果 明显 优 于 其 他 三 种 方法 ,从 MBF 参 数 图 中 


Phantom sSVD 


MBF 


MBV 


MTT 
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AY LAR EU TE aS LAL LID AL, EA 
可 以 发 现 ,本 文 方法 得 到 的 MBF,MBV 和 MTT 人 参数 图 
中 背景 区 域 和 心脏 区 域 条 形 伪 影 的 抑制 效果 显著 ;通过 
与 其 他 三 种 方法 对 比 , 本 文 方法 降低 噪声 同时 保持 精确 
定量 的 参数 图 估计 。 图 3 分 别 对 应 于 图 2 中 的 MBE， 
MBV 和 MTT 的 水 平 轮廓 图 (Profile)。 相 对 于 sSSVD， 
bSVD 和 Tikhonov 方 法 ,本 文 方法 所 得 的 轮廓 图 在 心 
室 , 缺 血 心 肌 和 健康 心肌 区 域 都 更 接近 于 体 模 参 数 图 的 
轮廓 图 ,继而 说 明 本 文 方法 能 够 实现 低 剂 量 CEMPI 血 
流动 力学 参数 图 的 优质 成 像 。 


Tikhonov 


图 2 不 同 去 卷 积 方法 得 到 XCAT 体 模 MBF、MBV、MTT 血 流动 力学 参数 图 


Fig.2 Hemodynamic parameter maps of the XCAT myocardial perfusion phantom by different methods. MBF in unit of mL-100 g`- min’, 


MBV in unit of mL/100 g, MTT in sec. 


为 了 定量 分 析 这 4 种 方法 的 差异 性 , 表 1 列 出 了 图 
2 中 各 参数 图 的 PSNR 和 rRMSE 测 量 值 。 由 表 1 可 以 看 
出 ,相对 于 sSVD,bSVD 和 Tikhonov 方 法 ,HDTV 方 法 


灌注 数据 仿真 猪 的 低 剂 量 心肌 灌注 数据 。 图 4 给 出 了 
猪 的 高 剂量 和 低 剂 量 仿真 心肌 灌注 图 像 。 
图 5 显示 了 不 同方 法 计算 的 心肌 血 流动 力学 参数 


可 以 获得 最 高 的 PSNR 测量 值 和 最 小 的 rRMSE 测 量 
值 ;尤其 从 MBF 参 数 图 的 定量 分 析 表 明 本 文 方法 获得 
更 高 的 PSNR 和 更 低 的 IRMSE。 
2.2 动物 数据 实验 

为 进一步 验证 本 文 方法 的 有 效 性 ,我 们 对 一 组 猪 的 
CT 心肌 灌注 数据 进行 了 测试 。 实 验 数据 由 GE 公司 
Discovery CT750 64 层 CT 设备 采集 , 射 源 到 扫描 物体 
和 探测 器 的 距离 分 别 为 538.52 mm 和 946.746 mm; 管 
电压 和 管 电流 分 别 为 120 kV, 100 mA, 层 厚 为 5 mm, 
CT 机 架 旋转 一 周 时 间 为 0.4 s; 旋 转角 在 [0,2x] 间 采样 
值 为 984 ,探测 器 个 数 888。 为 了 减少 辐射 剂量 避免 二 
次 扫描 ,根据 文献 [22-23 ] ,我 们 通过 获得 的 高 剂量 心肌 


图 。 由 图 5 可 以 看 出 ,本 文 提出 的 HDTV 正则 化 去 卷 积 
方法 抑制 了 噪声 和 伪 影 的 同时 保持 了 清晰 的 细节 。 从 
MBF 参 数 图 可 以 看 出 本 文 方法 估计 的 MBF 和 高 剂量 
得 到 的 MBF 具有 高 度 一 致 性 ,而 SSVD,bSVD 和 
Tikhonov 方 法 获得 的 MTT 参 数 图 严重 退化 。 基 于 奇异 
值 分 解 的 Tikhonov 方 法 获得 的 三 种 参数 图 在 噪声 较 少 
和 条 行 伪 影 抑制 方面 有 一 定 的 改善 ,但 是 通过 总 体 对 
比 ,本 文 方法 在 噪声 滤 除 和 边缘 细节 保持 两 方面 表现 更 
为 突出 。 表 2 所 示 为 图 5 中 各 参数 图 的 PSBNR 和 rRMSE 
测量 值 ,可 以 看 出 ,本 文 方法 较 其 他 方法 在 真实 动物 数 
据 实验 中 具有 更 好 的 噪声 抑制 性 能 。 
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| Fig.3 Horizontal profiles (186 th row, 250 th to 
336 th columns) of hemodynamic parameter 
maps in figure 2. A: Result from MBF; B: 
Result from MBV; C: Result from MTT. 
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表 1 对 应 图 2 中 各 参数 图 的 PSNR 和 rRMSE 
Tab.1 PSNR and rRMSE measured on the parameter maps in Fig.2 


MBF (mL: 100 g'+ min") MBV (mL/100 g) MTT (s) 
Method 
PSNR rRMSE PSNR rRMSE PSNR rRMSE 
sSVD 22.2456 0.4638 19.1955 0.6438 6.9944 4.7795 
bSVD 21.7904 0.6251 21.9737 0.4823 4.6393 6.3393 
Tikhonov 22.1086 0.5563 21.3343 0.5352 6.3917 5.5313 
HDTV 34.0013 0.087 32.5322: 0.1459 19.2882 1.2750 


图 4 高 剂量 和 低 剂量 仿真 猪 心肌 灌注 
CT 图 像 

Fig.4 Preclinical porcine high-dose 
myocardial perfusion image and 
low-dose image reconstructed from the 
simulated projection data. A: High- 
dose CT-MP image reconstructed by 
the FBP; B: Low-dose CT-MP image 
reconstructed by FBP. 


3 结论 验 。 通 过 与 SVD .bSVD 和 Tikhonov 去 卷 积 方法 实验 

针对 低 剂量 CT 心肌 灌注 去 卷 积 估计 动力 学 参数 结果 对 比 ,本 文 方法 实验 结果 表明 能 有 效 抑制 心肌 血 流 
不 准确 问题 ,本 文 提出 一 种 基于 HDTV 正则 化 的 去 卷 积 。 ”动力 学 参数 图 噪声 和 伪 影 , 且 本 文 方法 的 实验 结果 更 加 
方法 ,其 中 ,HDTYV 正则 化 利用 心肌 血管 空间 结构 相似 ”逼近 XCAT 真 值 参数 图 和 高 剂量 猪 参 数 圆 。 我 们 从 
性 和 心肌 血 流 信号 时 间 连 续 性 。 为 了 评价 本 文 方法 的 ” MBF,MBV 和 MTT 人 参数 图 可 以 清晰 判断 出 缺 血 心肌 的 
优异 性 ,我 们 采用 XCAT 仿 真 和 猪 心肌 灌注 数据 进行 实 。 ”位 置 和 大 小 ,这 对 于 临床 有 实际 的 应 用 意义 。 
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图 5 不 同 去 卷 积 方法 得 到 猪 的 MBF.JMBV、MTT 血 流动 力学 参数 图 


Fig.5 Hemodynamic parameter maps of pig myocardial perfusion image by different methods. MBF in unit of mL/100 g/min, 


MBYV in unit of mL/100 g, and MTT in sec. 


表 2 对 应 图 5 中 各 参数 图 的 PSNR 和 rRMSE 


Tab.2 PSNR and rRMSE measured on the parameter maps in Fig.5 


MBF (mL/100 g/min) MBV (mL/100 g) MTT (s) 
Method 
PSNR rRMSE PSNR rRMSE PSNR rRMSE 
sSVD 23.6417 0.4134 29.2716 0.3206 15.3749 0.1694 
bSVD 21.1776 0.6426 28.3514 0.4505 16.6818 0.2154 
Tikhonov 21.2826 0.6753 27.9775 0.4631 17.5036 0.3517 
HDTV 34.7795 0.1317 34.6883 0.1805 21.0989 0.0865 


由 于 算法 引入 惩罚 参数 ,来 平衡 保 真 项 和 正则 化 
项 ,本 文 未 对 4 进行 目 适应 优化 设计 ,而 是 通过 去 卷 积 过 
程 的 视觉 评价 ,因此 存在 一 定 的 缺陷 ;另外 ,本 文 仅 列 出 
了 三 种 商业 软件 广泛 使 用 去 卷 积 方法 作为 对 比 ,还 有 其 
他 重建 和 图 像 后 处 理 方法 需要 进一步 的 对 比试 验 ;还 有 
本 文 实验 基于 XCAT 体 模仿 真 和 动物 实验 ,名 要 实现 本 
文 方法 的 临床 应 用 ,还 需要 进行 更 为 广泛 的 临床 实验 。 


参考 文献 : 

[1] Becker A, Becker C. CT imaging of myocardial perfusion: possi- 
bilities and perspectives| J]. J Nucl Cardiol, 2013, 20(2): 289-96. 

[2] George RT, Jerosch-Herold M, Silva C, et al. Quantification of 
myocardial perfusion using dynamic 64-detector computed tomo- 
graphy[J]. Invest Radiol, 2007, 42(12): 815-22. 

[3] Mehra VC, Ambrose M, Valdiviezo-Schlomp C, et al. CT-based 


myocardial perfusion imaging-practical considerations: acquisition, 


image analysis, interpretation, and challenges [J]. J Cardiovasc 
Transl Res, 2011, 4(4): 437-48. 

[4] Brenner DJ, Hall EJ. Computed tomography-an increasing source of 
radiation exposure[J]. N Engl J Med, 2007, 357(22): 2277-84. 

[5] Speidel MA, Bateman CL, Tao Y, et al. Reduction of image noise in 
low tube current dynamic CT myocardial perfusion imaging using 
HYPR processing: a time-attenuation curve analysis[J]. Med Phys, 
2013, 40(1): 011904. 

[6] Ramirez-Giraldo JC, Yu L, Kantor B, et al. A strategy to decrease 

partial scan Reconstruction artifacts in myocardial perfusion CT: 

phantom and in vivo evaluation| J]]. Med Phys, 2012, 39(1): 214-23. 

[7] Tao Y, Chen GH, Hacker TA, et al. Low dose dynamic CT 

myocardial perfusion imaging using a statistical iterative Reco- 

nstruction method[J]. Med Phys, 2014, 41(7): 071914. 

[8] Biao ZY, Gong CF, Zeng D, et al. Motion Adaptive Sparsity Prior 


for Low-dose Dynamic CT Myocardial Perfusion Imaging[D]. The 
13th International Meeting on Fully Three-Dimensional, 2015. 


[9] Frindel C, Robini MC, Rousseau D. A 3-D spatio-temporal deco- 


201712.00732v1 


chinaXiv 


http://www.j-smu.com 


nvolution approach for Mr perfusion in the brain{J]. Med Image 
Anal, 2014, 18(1): 144-60. 

[10] Fang R, Zhang S. Chen T, et al. Robust low-dose CT perfusion 
deconvolution via tensor total-variation regularization [J]. IEEE 
Trans Med Imaging, 2015, Feb 20. [Epub ahead of print] 

[11] Ostergaard L, Weisskoff RM, Chesler D, et al. High resolution 
measurement of cerebral blood flow using intravascular tracer 
bolus passages. Part I: Mathematical approach and statistical 
analysis|J |. Magn Reson Med, 1996, 36(5): 715-25. 

[12] Wu O, Ostergaard L, Weisskoff RM, et al. Tracer arrival timing- 
insensitive technique for estimating flow in Mr perfusion-weighted 
imaging using singular value decomposition with a block-circulant 
deconvolution matrix J]. Magn Reson Med, 2003, 50(1): 164-74. 

[13] Zierler K. Indicator dilution methods for measuring blood flow, 
volume, and other properties of biological systems: a brief history 
and memoir|J]. Ann Biomed Eng, 2000, 28(8): 836-48. 

[14] Fieselmann A, Kowarschik M, Ganguly A, et al. Deconvolution- 
based CT and MR brain perfusion measurement: theoretical model 
revisited and practical implementation details [J]. Int J Biomed 
Imaging, 2011: 467563. 

[15] Fahmi R, Eck BL, Vembar M, et al. Dynamic CT myocardial 
perfusion imaging: detection of ischemia in a porcine model with 
FER verification| J]. Proc SPIE, 2014: 9038. 

[16] Fahmia R, Eck BL, Vembar M, et al. Dose reduction assessment in 


J South Med Univ, 2015, 35(11): 1579-1585 


ChinaXiv (ERAT 


+ 1585 - 


dynamic CT myocardial perfusion imaging in a porcine balloon- 
induced-ischemia model[J]. Proc SPIE Med Imaging, 2014: 9033. 

[17] Calamante F, Gadian GD, Connelly A. Delay and dispersion effects 
in dynamic susceptibility contrast mri: simulations using singular 
value decomposition[J |. Magn Reson Med, 2000, 44(3): 466-73. 

[18] Mouannes-Srour JJ, Shin W, Ansari SA, et al. Correction for 
arterial-tissue delay and dispersion in absolute quantitative cerebral 
perfusion DSC Mr imaging [J]. Mag Reson Med, 2012, 68(2): 
495-506. 

[19] Hestenes M. Multiplier and gradient methods [J]. J Optim Theo 

Appl, 1969, 4(1): 303-20. 

[20] Yang J, Zhang Y. Alternating direction algorithms for L1-problems 

in compressive sensing [J]. SIAM J Sci Comput, 2011, 33(1): 

250-78. 

[21] Segars WP, Sturgeon G, Mendonca S, et al. 4D XCAT phantom for 

multimodality imaging research [J]. Med Phys, 2010, 37(9): 

4902-15. 

[22] Ma JH, Liang ZR, Fan Y, et al. Variance analysis of x-ray CT 

sinograms in the presence of electronic noise background|J]. Med 

Phys, 2012, 39(7): 4051-65. 

[23] Ma J, Zhang H, Gao Y, et al. Iterative image Reconstruction for 


cerebral perfusion CT using a pre-contrast scan induced edge- 


preserving prior[J]. Phys Med Biol, 2012, 57(22): 7519-42. 


(编辑 : 吴 锦 雅 ) 


